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1. INTRODUCTION 

The use of PV as a conventional resource will be possible in the future. PV can be combined with 
other clean energy sources, e.g., fuel cells and batteries [1]. Power electronic converters, which have recently 
made significant progress, are characterized by their high efficiency, high reliability, and merging different 
renewable energies [2]. The majority of the studies are on multi-input DC-DC converters in the field of 
power electronics. Ahrabi et al [3], Furkan Akar et al [4], Banaei et al [5], Danyali et al. [6], Ghavidel et al. 
[7], Parham Mohseni et al [8], Neng Zhang et al [9] have introduced a new, highly-efficient multi-input DC- 
DC converter suitable for renewable energies. The application of power electronic converters is expanding in 
modern power systems [10]-[16] due to the integration of renewable energy. While controlling transmission 
power had problems in conventional power systems, transmission power is controllable in modern systems. 
Problems associated with real-time control of transmission power are major factors leading to blackout in 
traditional power systems, especially when the grid faces cascade outages [17]-[24]. Nowadays, the lack of a 
suitable controller makes the system approaches to collapse. However, the control of transmission power 
between two busbars is almost impossible in conventional power systems. By employing power electronic 
converters, the transmission power can be controlled, and immediate action becomes possible for the 
operator. 
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This paper focuses on the control issues of transmission power. The flexibility and high reliability in 

transferring power are the features of the proposed system. An advanced DC-DC converter, three-phase 
inverter, and the multivariable decentralized control technique are utilized for this purpose. 
There is little research on controlling DC-DC converters connected to a three-phase inverter to transfer the 
maximum power of renewable sources. A suitable multivariable control technique should be adopted because 
power electronic converters are multivariable systems. The control method must extract the maximum power 
of PV and simultaneously control active and reactive power. Khaki-Sedigh, and Moaveni [25] has delineated 
different multivariable control techniques, and this paper has employed a decentralized control technique to 
control the proposed system. Figure 1 shows the proposed system. 

The operation of modern power systems has faced new challenges upon the expanded integration of 
PVs. On the one hand, the maximum power of PV must be extracted at any moment to be economically 
justified. On the other hand, uncertainty in PV generation will cause problems, including balancing in power 
generation and consumption. Therefore, herein, a hybrid structure consisting of a battery and fuel cell is 
utilized to control the transferring power proportional to the demanded power while the maximum power of 
PV is extracted. A compromise can be made between uncertainty in the generation power of PV and the 
power required by the grid via the proposed hybrid structure; if the output power of PV is more than the 
demanded power, the proposed controller must charge the battery, and if power demanded by the network is 
more than the output power of PV, the control system must add a fuel cell and battery to the structure to 
transfer the set power to the network. 

The proposed system and controller pursue several objectives simultaneously: 

a. Extracting the maximum output power of PV at any moment 

b. Managing transmission power with the help of the battery and FC (without disturbing the maximum 
power extraction of PV) 

c. Regulating DC-link voltage 

The important point in this structure is that the transmission power does not affect the extraction of the 

maximum power of PV. 

First, the dynamic model of a hybrid DC-DC converter is achieved. Then, a suitable controller is 
designed. Subsequently, the phase-shift technique is applied to a three-phase inverter to transfer the set 
power. Finally, extensive simulation is prepared to validate the proper performance of the proposed system. 
Furthermore, fast fourier transform (FFT) analysis is conducted in each mode to present detailed results of 
the proposed system. 


Three-input Three- 
Hybrid DC- phase 


DC inverter 





AC Side 


Figure 1. General scheme of the proposed system 


2. SMALL-SIGNAL MODELING OF HYBRID DC/DC/AC CONVERTER 

Figure 2 displays the hybrid DC-DC converter that is the inverter input. Since the output voltage and 
current of the DC-DC converter are DC, the inverter and grid or the AC side are modeled only with simple 
resistance R_=Vo/Io. However, the designed controller is applied to the true system to verify the 
effectiveness of the controller. According to Figure 2, it is clear that a discontinuous current is imposed on 
PV, which disturbs the maximum power point tracking. A conventional buck-boost is allocated at the 
entrance of V2; consequently, Ipv=IL2xd2. Thus, this problem can be solved by controlling I2; in this case, the 
maximum power of PV is followed. The proposed system has three modes. 

In the first mode, PV and FC can produce the desired consumption power, so the battery must be 
bypassed. In this mode, it is assumed that the requested power by grid is more than the maximum power of 
PV. In the second mode, the total generation power of PV and FC cannot supply the output, so the battery is 
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discharged. In the third mode, it is assumed that the whole generation power of PV and FC is larger than the 
required consumption power in the output stage, so the surplus energy is saved in the battery. A separate 
controller must be designed in each mode because the state-space model of the system is different. All three 
modes cover the reality in the grid. Therefore, in first step, the controller must identify the mode of operation 
and then in the second step, act such that the intended control objectives be achieved. 





Figure 2. The Hybrid DC-DC Converter 


2.1. First operation mode (Supplying the load without the contribution of the battery) 

To bypass the battery in this mode, S3 and S4 must be turned off and turned on, respectively. 
Consequently, d4 (duty cycle of S4) =1 and d; (duty cycle of S3) =0; thus, d4 and dz as control signals are 
missed. Only dı (duty cycle of S1) and d2 (duty cycle of S2) are utilized to extract the maximum power of PV 
and regulate the output voltage. Furthermore, (1) shows the state-space model of systems. 





diy, dV, l , , , 
1 a = —ryip, + (V1 + Vce)di + (V1 — Vo) (4 — dy), C—- = (ip2 — ip1)(di — d2) — it1d2 + ip2(1-d2) 
dilz ! dVo . V 
L2 = = —r2İ12+V2d2 — Ve(1-d2), Co w i11 (1-d1) — E (1) 
L 


The (1) is driven from the principle of the voltage balance and current balance in the inductor and capacitor 
in the steady-state, respectively. Moreover, (1) reveals that reaching the transfer function matrix is impossible 
because the duty cycle of switches is directly related to state variables. The linearization technique at the 
operating point is a well-known technique for solving this problem. In this method, state variables, duty 
cycles, and inputs are divided into two component DC values (X,V,D) and perturbation 


(%,¥,d): 
x=X+X%v=V4+i0,d=D+d (2) 
If it is assumed that the perturbation is small and does not significantly vary during the switching period 


(x << X ,¥<<V.d<<D), by substituting (2) into (1) and neglecting the second terms, small-signal models are 


represented in the matrix form as (3): 


x 
$=CX+Du (3) 
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where X, U, and ÿ are state variables' vector, control variables' vector, and the output system, respectively. 
Therefore, the matrix form of the small-signal model for the first mode is obtained as (4): 














dip Dı -1+D, YOUVE) 
dt Lı Ly Ly Lı 7 
dij2 022 -1+D2 0 Ti 0 V2tVc}__ Na 
dv. | ~ i y Tad P= g pT 
c “Di D3 00 Ve tui t2 || do 0001! | Ve 
dt C C To C C Vo 
dVo 1-D1 = ay 
Co RLCọo Co 


2.2. Second operation mode (Supplying the load with discharging the battery) 

In this mode, the maximum power of PV and FC must be extracted, and the battery is connected to 
the system to regulate the output voltage. S4 must be kept turned on to discharge the battery. Moreover, the 
power flow to the battery is controlled by turning on and turning off S3. Furthermore, (5) expresses the state- 
space model of systems. 


di, dv 
Lı —— =-r1i11 + (V1 +Vc+Vg)d; + (Vi +Vc)(d; — d3) + (Vy — Vo)(1-d,), oars oe + (irz — ini) (di — d2)+iL2(1-d1) 


di, 
L2 a r2iL2 + (V2+Vp)d3t+V2(d2 — d3)-Vc(1-d2z), Co e =i, (1- = (5) 


Similar to the first mode, the dynamic model is obtained as (6): 





ditt Dy -1+D1 Vo+Vc 0 VB 
dt 7 Lı ly y Ly Ly is 
diy.2 0 72 -1+D2 o ‘L1 025c% || d, 1000 m 
L2 L l L2 L x l 

dt 2 L2 2+. 2 2 11a | y=|0100|| 27]; D=0 (6) 
Ef 00 || 0 jla 0001] | <° 
dt C Ves C Cc ds Vo 
do 1-D af -I 
ETa = 00 — 00 

RCo Co 


2.3. Third operation mode (Supplying the load with charging the battery) 
The charge path is established by turning off the S4, and the current of the battery is controlled by 
turning on and off the S3. Moreover, (7) shows the state-space model of systems. 


d 
di No i s 
L2 = =-r2i12 + V2d3+(V2 — Vg)(d2 — d3)-Vc(1-d2), Co dt Ut Oa) Ri n 


di : dV. ' : : 
L = =Tyipy + (Vi +Ve)d3 + (Vi +Ve — Ve)(di — d3) + (V1 — Vo) (4-d1), C z =-inide + Giz — in) (di — d2)+ir2(1-d1) 
dt 


Similar to the previous mode, the dynamic model is obtained as (8): 





ditt -r1 n Da -14Di Vo+Vc-Vp 0 VB 
dt Ly Ly Ly Lı Lı 
dip2 o 2 ri+D2 0 tha o W2tVc-Ve Ve d, 1000 tha 
dt Lo L 12 L Lo ||5 112 
P AN aea as * 1/d,|.y=]0100]}~“ |; D=0 (8) 
a | [00 |] 0 : 0001) |° 
dt C V C C d3 Vo 
dVo 1-D =j -I 
a 11 99 00 
Co RLCo Co 


3. DESIGNING THE CLOSED-LOOP CONTROLLER AND EXAMINING ITS PERFORMANCE 
In the first mode, the main control objectives were specified. IL2 is controlled to PV working at 
maximum power, and the output capacitor voltage is controlled to adjust the output voltage. As explained 
before, in the first mode, the duty cycle of S3 and S4 is missed as control signals, and only dı and dz are 
available as control signals due to bypassing the battery. Now, one must control state variables (IL2, Vo) with 
dı or d2. In multivariable control, control signals for each output (state variables) are selected by the relative 
gain array (RGA) matrix. The RGA is a widely-used classical method for determining the best input-output 
paring for the multivariable process control system. In its general form, the RGA matrix is defined as (9): 
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The RGA matrix is calculated in the steady state as: 


; Zo) =L39 

A(G(0))=G(0)*670) = (“F355 539) (10) 
2.1365 —3.2139 2.0774 

NG(O)-GO).*6" ={ 0 3.3333 -23333 (11) 
—1.1365 0.8806 = 1.2559 
2.2656 —1.4350 0.1694 

Neonato =È 0 1.8224 -08224 (12) 
—1.2656 0.6126 1.6530 


(10), (11), and (12) are the RGA matrix for the first, second, and third modes, respectively. 

The RGA matrix shows the dependency of the control signal to the output. In (10)-(12), diagonal 
elements have significant values compared to other elements, so it can be concluded that, in (10), dı and d2 
have the most impact on Ir2 and Vo, respectively. However, being positive has a priority to the value in the 
RGA matrix. Therefore, pairing the input-output set is defined as follow: 


First mode: [,, >d; Vo>d, 
Second and Third modes: J,, >d, I >d, Vo >d, 


Figures 3 and 4 illustrate the closed-loop configuration in the first, second, and third mode. The 
most important part of the multivariable system is system coupling. In these schemes, the CP matrix is placed 
to decouple the system, and PI controllers undertake reference tracking. There are several options for 
designing the CP matrix, but in general form, CP depends on systems and the control objective. The CP is 
chosen as follows for decoupling the system in the steady state: 


cP=[G(0)]}" (13) 





Figure 4. The closed-loop scheme for second and third modes 
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The closed-loop performance is simulated with the help of MATLAB/SIMULINK. Figure 5 shows 
the closed-loop performance in the first mode. In Figure 5, U = [4 0] is applied to the system in the second 
moment. It is expected that Yı and Y2 should have value due to the inherent coupling in multi-input systems. 
The CP matrix makes the system decouple in the steady state, and the PI controller causes reference tracking 
to occur. The accurate performance of the controller is evident because it managed to eliminate the intrinsic 
coupling in the system, and reference tracking happened simultaneously. All these features are achieved by a 
simple CP matrix and diagonal PI controller. Figures 6 and 7 demonstrate closed-loop performance in the 
second and third modes. Similar to the first mode, good performance has been achieved. 
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Figure 5. Closed-loop performance in the first mode Figure 6. Closed-loop performance in the second 
mode 
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Figure 7. Closed-loop performance in the third mode 


4. SIMULATION RESULTS 

The closed-loop controller was designed in previous sections. In this section, the hybrid DC-DC 
converter is connected to a three-phase inverter to transfer the power of PV, FC, and battery to the grid. All 
modes are simulated by MATLAB/Simulink. Figures 8 and 9 display PV and FC characteristics, 
respectively. The ready block of PEMFC is a preset model for fuel cells in MATLAB/Simulink and has been 
employed in this paper. 

The simulation consists of three modes. In the first mode, the sun irradiation of PV is S=500W/m7, 
and under this condition, the maximum generation power of PV is about 10800 W, and it occurs at 300 V. 
Moreover, in this mode, the maximum power of FC equals 8325 W, and it happens at around 225 A. The 
battery is bypassed in this mode, and the total generation and transmission powers are assumed to be equal 
(10800+8325=19125 W). 

In the second mode, the sun irradiation of PV is S=600 W/m7?, and the maximum power under this 
irradiation is 12950 W and takes place at 300 V. FC has 6000 W in the second mode, and under this power, 
FC can give nearly 133 A. The summation of powers from both sources equals 12950+6000=18950 W; 
therefore, in comparison to the first mode, the total generation power has decreased. Consequently, the 
battery must be added to the system to resolve the lack of power and converter losses. 

In the third mode, the irradiation of PV is S=750 W/m’, the maximum power of PV reaches about 
16102 W, and the voltage of PV is close to 285 V. The power of FC is approximately 4284 W, and under this 
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power, FC can give 86 A. The total generation power from both sources equals 16102+4284=20386 W, 
which is larger than the demanded power for transferring. As a result, the battery must be charged to save 
surplus energy in the third mode. 

The state variables of the hybrid DC-DC converter are controllable. Hence, I1, IL2, and Vo must be 
controlled to extract the maximum power of resources besides regulating the DC-link voltage. The output 
voltage of the inverter depends on the DC-link voltage. Therefore, the inverter does not need to be controlled 
in a conventional way such as sinusoidal pulse width modulation. Furthermore, the output voltage of the 
inverter can be controlled by changing the DC-link voltage (Van=2Vac/n). The only parameter that must be 
controlled in the inverter is the phase difference. In all three modes, the inverter must be able to produce a 
+90° phase difference (6 = 90° ); in this case, maximum active power is transferred from the DC side to the 
AC side. 
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Figure 8. PV characteristics Figure 9. Fuel cell characteristics 


A 6-s simulation with three power management strategies is provided to evaluate the proposed 
control system in each operation mode. Figure 10 and Figure 11 depict the duty cycles of S1, S2, S3, and S4, 
respectively, and show the control signals that vary to extract the maximum power of resources. Figure 12 
represents the hybrid DC-DC converter outputs. The most important part of the control system is the duty 
cycles of the converter switches, which have been shown in Figure 10 and Figure 11 because the transient- 
state of the converter depends on the variation of duty cycles. 
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Figure 10. Duty cycles of S1 and S2 Figure 11. Duty cycles of S3 and S4 


Based on Figure 12, it is confirmed that FC works at maximum power due to I,;=225A in the first 
mode, ILı=133A in the second mode, and ILı=86A in the third mode. Since Iu; is the exact Ipc, as expected, 
the FC operates at maximum power. The current of first (IL1) and second (112) sources, output voltage of DC- 
DC converter (Vo), and the current of the battery (Ibattery) have been shown in Figure 12. Proper transient- 
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state, suitable overshoot and undershoot of DC-DC converter state-variables, and reference tracking can be 
clearly seen in Figure 12, which is the result of the appropriate operation of the control system. 

Figure 13 depicts the output voltage of the inverter. Figure 14 verifies that PV operates at maximum 
power. To prove that the inverter managed to produce a +90° phase difference in each mode, the output 
voltage of the inverter must be examined by FFT analysis. In all the modes, the voltage amplitude of the grid 
is assumed to be 100 V with a 60Hz frequency. Since all the calculations are based on amplitude, a 3/2 factor 


appears in basic equations. 
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Figure 12. Output of the hybrid DC-DC converter 
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Figure 13. Output voltage of the inverter Figure 14. Output voltage and current of the PV 


First-mode simulation 0 <t < 2 ; Figure 15 represents the output voltage of the inverter. As expected, 
the fundamental harmonic has a +90° phase difference, and it causes the maximum power of PV and FC to be 
transferred to the AC side or gird side. A simple phase shift can produce the desired phase difference in the 
inverter. Figure 16 shows the FFT analysis of the current of the inverter (Ian). 

FFT analysis in Figure 15 (a) and Figure 15 (b) confirms that the inverter produces a +90° phase difference, 
according to: 


3ViV2.. V1=300.3V2=100XL=2.64Q 
P.2(39) = 2 Xi: sin 6 ——— Peo = 17069W 





On the other hand, FFT analysis gives the phase differences between current and voltage: 
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VI V=300.32+90°I=120218.4° 
Å= a 


Sin v(30) — 2 


3 
Sin v(3q) => X (300.3) x (120)290° — (180 + 18.4)° 


Sinv(39) = — 17063 — j51290.5 

Pi2 confirms Siny, Which means that the inverter transfers 17063 W of active power to the grid. Only a 6W 
difference between Siny and P12 is originated from numerical errors, which is acceptable in this range of 
power. 
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component amplitude component amplitude 


Second-mode simulation 2 < t < 4; Figure 17 (a) and Figure 17 (b) presents an FFT analysis of 
Van. It verifies that the inverter produces a +90° phase difference. Similar to the first mode, Figure 18 (a) and 
Figure 18 (b) shows the FFT analysis of the current of inverter (Ian). 


3V1V>_. V1=325.32+909V2=100209 XL=2.64 
Pi2(39) = os sind — a = 18490.4 


Sinv can be calculated from the FFT analysis in the second mode according to Figure 17 and Figure 18: 
Van=325..32+90°ian=129.1217.1° 
Sin v(39) 


=r C—O 
ore == x 325.3 x 129.12 — 107.1° 


Sinv(3q) = — 18523 — j60209.56 
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Sinv and Pı2 confirm each other well. A 32W difference is originated from numerical errors, which is 
acceptable in this range of power. In all the equations, V;=Van (the inverter voltage) and V2=V gira. 
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Figure 17. FFT analysis of Van in the second mode, 
(a) inverter output voltage; (b) fundamental 
component amplitude 
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Figure 18. FFT analysis of Ian in the second mode, 
(a) inverter output current; (b) fundamental 
component amplitude 


Third-mode simulation 4<t <6 ; Figure 19 indicates the power of the battery in three modes. Similar to the 


first and second modes Figures 20 (a), (b) and Figures 21 (a), (b) depict the FFT analysis in the third mode. 
Figure 20 proves that the inverter creates a +90° phase difference. 
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Figure 19. Power of the battery in all the modes 
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Figure 20. FFT analysis of Van in the third mode, 
(a) inverter output voltage; (b) fundamental 
component amplitude 


Figure 21. FFT analysis of Ign in the third mode, 
(a) inverter output current; (b) fundamental 
component amplitude 


5. CONCLUSION 

A new multivariable control technique was employed here to transfer the intended power to the grid 
side and extract the maximum power of PV simultaneously. The performance of the controller was examined 
under different scenarios. Simultaneous active power control and DC-link voltage regulation were possible 
for the proposed system. 

In the proposed system, the three-phase inverter worked with a 50Hz switching frequency, due to 
which the leakage current of PV was no longer a problem. Power management and online operation can be 
performed more reliably in micro-grids and modern power systems with the proposed system. Still, more 
research on control and various power management strategies 1s required. 
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